home *** CD-ROM | disk | FTP | other *** search
/ Shareware Grab Bag / Shareware Grab Bag.iso / 004 / mathstat.arc / S13.BAS (.txt) < prev    next >
Encoding:
GW-BASIC  |  1982-04-22  |  1.6 KB  |  83 lines
  1. 10   PRINT "MULTIPLE LINEAR REGRESSION"
  2. 20   PRINT 
  3. 30   DIM X(9),S(9),T(9),A(9,10)
  4. 40   PRINT "NUMBER OF KNOWN POINTS";
  5. 50   INPUT N
  6. 60   PRINT "NUMBER OF INDEPENDENT VARIABLES";
  7. 70   INPUT V
  8. 80  X(1) = 1
  9. 90   FOR I = 1 TO N
  10. 100   PRINT "POINT ";I
  11. 110   FOR J = 1 TO V
  12. 120   PRINT "  VARIABLE ";J;
  13. 130   INPUT X(J + 1)
  14. 140   NEXT J
  15. 150   PRINT "  DEPENDENT VARIABLE";
  16. 160   INPUT X(V + 2)
  17. 170   FOR K = 1 TO V + 1
  18. 180   FOR L = 1 TO V + 2
  19. 190  A (K,L) = A (K,L) + X (K) * X (L)
  20. 200  S (K) = A (K, V+2)
  21. 210  NEXT L
  22. 220  NEXT K
  23. 230  S (V + 2) = S (V+2) + X (V+2) ^2
  24. 240  NEXT I
  25. 250  FOR I=2  TO V+1
  26. 260  T (I) = A (1,I)
  27. 270  NEXT I
  28. 280  FOR I=1  TO V+1
  29. 290  J = I
  30. 300  IF A (J,I) <> 0  THEN 340
  31. 305  J = J + 1
  32. 310  IF J < = V + 1  THEN 300
  33. 320  PRINT "NO UNIQUE SOLUTION"
  34. 330  GOTO  810
  35. 340  FOR K = 1  TO V + 2
  36. 350  B = A (I,K)
  37. 360  A (I,K) = A (J,K)
  38. 370  A (J,K) = B
  39. 380  NEXT K
  40. 390  Z = 1 / A (I,I)
  41. 400  FOR K=1  TO V+2
  42. 410  A (I,K) = Z * A (I,K)
  43. 420  NEXT K
  44. 430  FOR J=1  TO V+1
  45. 440  IF J=I  THEN 490
  46. 450  Z =-A (J,I)
  47. 460  FOR K=1 TO V+2
  48. 470  A (J,K) = A (J,K) + Z * A (I,K)
  49. 480  NEXT K
  50. 490  NEXT J
  51. 500  NEXT I
  52. 510  PRINT
  53. 520  PRINT "EQUATION COEFFICIENTS:"
  54. 525  PRINT "     CONSTANT:"; A (1,V+2)
  55. 530  FOR I=2  TO V+1
  56. 540  PRINT "VARIABLE("; I-1;"): ";A (I,V+2)
  57. 550  NEXT I
  58. 560  P = 0
  59. 570  FOR I =2  TO V+1
  60. 580  P = P +A (I,V+2) * (S(I) - T (I) * S (1)/N)
  61. 590  NEXT I
  62. 600  R=S (V+2) - S (1) ^2/N
  63. 610  Z = R-P
  64. 620  L = N - V -1
  65. 640  PRINT
  66. 650  I = P/R
  67. 660  PRINT "COEFFICIENT OF DETERMINATION (R^2) ="; I
  68. 670  PRINT "COEFFICIENT OF MULTIPLE CORRELATION =; SQR (I)
  69. 680  PRINT "STANDARD ERROR OF ESTIMATE ="; SQR (ABS(Z/L))
  70. 690  PRINT
  71. 700  PRINT "INTERPOLATION: (ENTER 0 TO END PROGRAM)"
  72. 710  P = A (1, V+2)
  73. 720  FOR J=1  TO V
  74. 730  PRINT "VARIABLE" ;J;
  75. 740  INPUT X
  76. 750  IF X=0  THEN 810
  77. 760  P = P+A (J+1, V+2) *X
  78. 770  NEXT J
  79. 780  PRINT "DEPENDENT VARIABLE="; P
  80. 790  PRINT
  81. 800  GOTO 710
  82. 810  RUN "STAT"
  83.